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HATIOHAL ADVISQHY' CCtMTETSE ECS AERQHAtJTICS 


EESEAECH MmCORANISM 



A STOUr CE SEEN TEWPERATURES OP COIIIC3AL BODIES 
IN SOPERSOHIC PEJX5EP 

By WllTser B. Hustoa, Calvin N. Warfield^ and Anna Z. Stone 


SUMMARY 


A comparison. Is made "between tiie time history of skin tompeirature 
measured on the nose of a "V-a and the bemperature computed using Eber's 
experimenteO. relation for haat- transfer coefficients for conical "bodies 
under supersonic conditions. The agreement obtained is felt to Justify 
the tise of B"ber's relation in the calculation of skin tenperatures under 
flight conditions. A general method developed for making such skin- 
tempera ttire calculations is used to oonrpute the variation of skin 
tenperature with time for a wide range of values of the pertinent 
parameters. The results show that "by proper selection of the "basic 
petrameters the Increase of akin temperatur*e during a limited time of 
flight can "be held to structurally permlssl"ble values. Methods a 3 ?e 
given for taking into account, when necessary, the effects of solar 
heating and the radiation exchanged "between the skin and atmosphere. 

Time histories of skin temperature are conputed for representative 
flight plans of a typical supersonic airplane. 


INTRODUCTION 


Because of the high stagnation and boundary “layer tenperatures 
associated with supersonic speeds, the effects of aerodynamic heating 
have been of considerable concern to tte designers of supersonic 
airborne structures. Values of akin temperature under steady- state 
conditions were calculated in reference 1. This report showed the 
necessity of allowing for the variable spooiflo heat of air in stagnation 
temperature calculations and by allowing for reidiation showed that 
skin temperatuireB could be lower than boundary -layer temperatures. The 
skin-temperature values calculated were high, however, from a stiructural 
or from an. operational steundpolnt, and the results were based on an 
extension of subsonic heat-transfer coefficients to a supersonic wedge- 
shaped airfoil. For this extension no esperimentai verification is 
yet available. 

Experimental Information on aerodynamic heating effects under 
supersonic conditions has been fragmentary. During the course of a 
series of T-2 flights at White Sands, New Mexico, the Naval Research 
Laboratory has, however, made a number of moaaurements of the variation 
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of skin ten^ieratiire with time . The data obtained on March 7, 19^7 on 
the conical warhead of missile number 21 ai:e of particular interest 
since these data can he ccm^ared with the results cf ceLLculatlons based 
on an experimental relation for supersonic heat-transfer coefficient 
for conical bodies given In reference 2. 

Since the resvilts of refei?enoe 2 are In such a form that they 
could be used in general calculations of the variation of skin tempera- 
ture with time, such a comparison Is of particular interest, especially 
if it Indicates that the results of reference 2 can bo extended to 
flight conditions. 

It Is the purpose of the present report to show the order of 
agreement between the time history of skin temperatures measured on 
the V-2 and the time history of temperatures calculated by use of the 
results of refez*ence 2. A general method of making such calculations 
Is given which makes It possible to calculate the time hlstozy of skin 
tenperature for a wide variety of structural euid fll^t-path parameters. 
The eq[uations given are sttfficlently general that the effects of solar 
radiation can be Investigated as well as the radiation received from 
space and the ambient atmosphere . 


SYMBOLS 


A area, ft^ 

C solar constant, 0.1192 Btu/ft^ sec 

c specific heat of skin material, Btu/lb ^ 

Op specific heat of air at constant piressuro, Btu/lb °F 

Fg skin factor (equal to ctw), Btu/ft^ °F 
g acceleration due to graviiy, 32.17^0 ft/sec^ 
h heat-transfer coefficient, Btu/ft^ sec °F 

h^ enthalpy per unit mass of air corresponding to total temperature, 
Btuyib 

hft enthalpy ner unit mass of air corresponding to ambient temperature, 

BtuAlj' 

H altitude, ft 

J mechanical equivalent of heat, 778.27 ft-lb/Btu 
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K teii^eratia:^ recovery factor 

k thermal oonluctlvlty of air, Btu/ft^ sec (^/ft) 

1 oharacterlstic length, ft 

Nu Nusselt number (hl/k) 

Pr Prandtl number (cp(ig/lc) 

p pressure of ambient atmosphere, Ib/ft^ 

Pq pressure of standard atmosphere at sea level (2116.229 Ib/ft^) 

Q quantity of heat, Btu 

R© Reynolds nimiber (VIp/p.) 

ambient atmospheric temperature, °F abs. 

% boundary -layer temperature, abs. 

Tg equilibrium skin tempera tiire , °F abs. 

Tjn mean skin tenqjeratiare over time interval At, °F abs. 

Tq temperature of standeqrd atmosphere at sea level ( 518 . 4 ° P abs.) 

Tp potential ten^erature , a fictitious temperattnre defined in 
equation (B2), °P abs. 

Tg skin temperature, °P abs. 

total or stagnation tenqeraturs , abs. 

t time, sec 

V velocity, ft/sec 

¥ weight, lb 

w specific weight of skin material, lb/ft3 
3 total apex angle of cone, radians 
7 adiabatic exponent 

5 sky radiation factor (See appendix A.) 

€ emlssivity 

P density of atmosphere, slugs/ft^ 
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Po 

a 

T 

0 


density of the standard atmosphere at sea IstsI (O. 002378 slTig/ft^) 
Stefan-Boltzmann radiation, constant (^.47594 x 10”^ 


Btu 


V 


ft® sec ahs 


■) 


sldn thiolmess, ft (or in.) 
viscosity of air, Ih-sec/ft® 


visooBlty of air at sea-level conditions (3 *7250 X lO”*^ Ih-sec/ft®) 
thermal lag consteint, sec 


BASIS FOR SKEN-^MEEEATUES CALCDLAa?IOHS 


Heat-halanoe equation .- Given the initial conditions and the 
specified fU^t path, the skin teniperature can he evalxxated as a 
function of time hy means of the following differential eq.uation which 
is derived in appendix A: 

^ ^ h(TR - Ts) + €o(&Fa^ - + gC 

dt o T w 


Equation (l) teOces into aocotmt the heat gained or lost hy thermal 
conductivity in the boundary layer, the exchange of radiant energy 
between the skin and the surrounding atmosphere, and the heat gained 
from the sun. The physical properties of the skin material are taken 
into account in the product or w which in the rest of this report will 
he called the skin factor Fg. The equation applies to a thin skin 
under the assuB^tlon that there is no conduction or radiation of 
heat to other parts of the structure . 

Eqxiatlon (1) is a nonlli»ar differential eqiiatlon for which 
the explicit Integral la not usefta, hut which can he integrated hy 
numerical nethods. Using average values of the various quantities 
over successive time intervals At, the relation between skin 
ten^rature and time can he expressed in the following equation; 

TSg " ’fel * - %) + « ^ - %'*]! ^ (2) 

The subscript 1 la used to denote the start of the Interval and 
subscript 2, the end of the Interval. The bar is used to denote 
average values in the time Interval At. 



mck BM Bo. L7J21 


5 


In. order to inake use of eq.uations (l) and (2), it is necessary 
to know the veO-xie of heat-toansfer coefficient h and boundary- 
layer tengperatxrre 3^. Methods for determining these factors, based 
on experiments at supersonic speeds, have been given in reference 2, 
and are utilized in the present report. 

Heat- transfer coefficient .- The value of h to bo used in 
equations (l) or (2) may be dstermlned from reference 2 in which 
it was shown that heat-transfer data could be correlated on the basis 
of nondimensional parameters and a shape factor by means of the equation 

Bu = (O.OOTl + 0.0154 (3) 


This eqxjation correlates wind-tunnel experiments on a series 
of conical models with apex euagle P of 10° to 120° covering a 
Mach number range of 1.2 to 3*1 Beynolds number range of 

2 X 10^ to 2 X 10^. In conputing the value of Bu and Be, the 
conductivity and coefficient of viscosity were based on the boundary- 
layer tenperature, the density lused was the density in the free stream 
ahead of the model, and the characte3^stic length was the total length 
of the conioal, model as meastzred along the surface (length of generatrix]^ 
and the value of heat-transfer coefficient was the average value for 
the entire surface. The Prandtl number was taken as constant and ignored 
in determining equation (3), and only e::qperimental data for which the 
shock wave was attache^ to the model were included. 

In the extension of equation (3) to the calculation of akin 
temperatuies in fUgbt, it is necessary to select a tenperature 
which characterizes conditions within the bou n d ar y layer, and also 
a characteristic length, in oirder to evaluate h. Ihe experimental 
data of reference 2 were obtained under transient conditions which 
differed from equilibrium by, at most, 36 ° S'* The characteristic 
tenperattire used as appropriate to such conditions is the boundary - 
layer temperature or stagnation teuperature suitably corrected by the 
recovery factor K. Under fli^t conditions with a finite thickness 
of structviral material, the skin temperature can be as much as 
several thousand degrees less than the boundary -layer tenperattire Tg. 

The skin temperature will modify the teuperature of the air In the 
Imnediate vicinity of Ihe skin and consequently the hsat -transfer 
capaoi-iy of the bovindary layer. There is, therefore, a question as to 
whether the value of h should be based on 3^, on or perhaps 

on some intermediate temperature . 

For skin teuperature calculations in which conduction of heat 
ai ong the thin skin may be neglected, as in equation (l), the local 
value of heat- transfer coefficient at any particular point is needed. 
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In the extension of eq^uation (3) to flight conditions It has heen 
somewhat arhltrarlly assumed that eqxiatlon ( 3 ) can he used to obtain a 
local value of h if the characteristic length is taJton as the 
distance of the particular point from the vertex of the conical body, 
measured along the surface . The selection of the proper characteristic 
length is aJLso a part of the q.uestion of the rangd of Reynolds number 
over which equation (3) may be considered valid. In the calctilatlona 
of the present report, equation ( 3 ) has been used regardless of the 
value of Re . In the absence of other experimental supersonic heat- 
transfer data, the validity of the assunptions about a cheLracteristlc 
temperature, chaia.cteristic length, and Reynolds number may be tested 
by coniparleon with the results of flight tests. 

. For convenieiace in oea.ctilation of the value of the heat'-transfer 
coefficient, equation ( 3 ) laay be rewritten in the following form. 



In ihe evalmtlon of equation (^), the values of jj, and k 
corresponding to values of IJg less them 2hOO^ F abs. may be based 
on table 3 ot reference 3* For values of greater than 2400° F abs., 

approximate values for ji and k may be determined using the assumption 
that the Prandtl number CpUg/k is equal to O .65 at high ten^ratures 
and that 


i - 



(5) 


Equation (5) ie a tentative approximation based on extrapolation of the 
data in reference 3 and should bo used only in -the absence of experimental 
data on the viscosity of air at high temoperatures. 

Value of p/pQ for pressuare altitudes frcan sea level to 65^000 feet 
may be taken from tables of the NACA standard atmosphere (reference 4-); 
values of p/Pq for altitudes above 65,000 foot corresponding to the 
tentative standard atmosphere may be found in table V of reference 5 « 

For einy peu:ticular configuration (p and I specified) when the rolatlcn- 
ship between pressure, temperattu^ , and altitude is known (from 
measurements or from tables of the standard atmosphere), tbe value of 
the heat “transfer coefficient h can be computed when sQ.tltude and 
velocity (or Mach number) are specified. 

Boundary- layer tenperature . - The boundary-layer temperature 
may be evaluated from the equation 
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Tg = Ta + - 

Values of tih© temperature recovery factor K ^ for conical todies In 
supersonic flow are given in figure 1. | 

The stagnation ten 5 >erat'ure T>p correspoi^ding to any specified 
values of velocity and ambient temperature is -given ty the following 
relation for adiabatic ccsapresslble flow at constant total energy 



( 6 ) 


r^T 

Jta 



V dv 


0 


(7) 


If the range of temperature change Is such that Cp may be 
considered constant, eq,uatlon (7) can be Integrated directly, and the 
result of the integration can be expressed as 


% 



( 8 ) 


Although this expression Is a conrenlent one, values of Tiji given hy 

eq.uation (8) are too large at hl^ values of M. Eg^uatlon (7) iney, however, 
he conveniently evaluated hy use of table 1 of reference 3 which gives 
values of the function 


/ f T 

c dT 
400 ^ 

for the enthalpy or total heat of air over a range of tenperature T 
from 300° F abs. to 6500° F abs* In iising this table, the stagnation 
temperature may be obtained firom the relation. 


h^ 


= lu + 

2gJ 


(9) 


Values of T^i from eq.uation (9) ar» conpared with those given by 
equation (8) with 7 = 1.4 in table I. At a Mach number of o and 
Ta = 392 . 4 ° F abs., the stagnation temperature is 732° F less th an the 
value given by equation (8) with 7 = 1.4. 

Skin factor .- Since the thickness, specific heat, and specific 
weight of the skin material enter the differential equation for akin 
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ten^eratvire as the product of the three quantities, the resxilts cf an 
integration of equation (2) for a particular value of Fg will apply 

to any structural material of hlg^ thermal conductivity and of stifficient 
thlclmess to give the specified value of Fg. Values of c anfl v for 
several commonly used aircraft structural materials are given in 
table U. The valvies of 0 for vrou^t iron and stainless steel are 
average values for which the tes^rature rangs is not specifically 
given in reference 6. In view of the large tengjerature variation of the 
specific heat of iron, allowance may scanetlmes need to he made for the 
change of Fg with tes^erature when calculations are specifically 
concerned with iron or steel. (See 'reference 7.) Values of Fg corre- 
sponding to various thicknesses of aircraft structural materials (and, 
for iron and steel to various ten^ratures) are shown la figure 2. 


« CALCULATION CEF SFIN T^MFEMTURE! 


The application of equations (1) or (2) to the determination of 
skin temperature requires that certain basic parameters and initial 
conditions be specified. These basic parameters which characterize 
the flight plan and the structure are, 

M Mach number 

E altitude 

# 

I length 

p cone angle 

e emlssivity 
Fg skin factor 

In addition, when skin tesgterature is to be calculated as a function of 
time, the initial skin temperature must 'be specified. 

Time-history calculaticxtiB . - The variation of skin temperature 
with time while an airplane or missile fo3J.ows a particular flight 
path in which either ^eed or altitude or both vary is often of considerable 
Interest. Calciilatlons of such a time history are eaelly performed, 
using equation (2) in which average valvies of h, T^^ 5, and T;;^ are 

used. For preliminary design purposes, to determine the range of skin 
temperatures expected eind to determine the relative Inportanoe of the 
various basic parameters in limiting the tempera tvire to structurally 
possible values during a flight of limited duration, it is convenient to 
use a hypothetical flight plan correapondlng to the step function 
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freq^uently xised. In the anaJj-sls of transient phenomena. In this 
hypothetical fli^t pla n it is asstanad that the altitude is constant, 
that the value of M is reached Instantaneously at t = 0 , and that 
the initial skin tenperature is eg.ual to the temperature of the amblenb 
atmosphere . Tlme-hlst<n:y calculations usin^ this f li^t plan are 
particularly slsple; an outline of the method and an example are given 
in appendix B. 

Eq.uill'brium skin temperat^lre . - The oqtaollihrlxm skin temperature Tg 
is the skin teaperatirre which would he reached after a sufficient lapse 
of time vmder steady fli^t conditions at constant altitude . The value 
of Te is independent of F3 and may he approximately eq[ual to Tg, or 
very much less, depending on the value of h and the importance of the 
radiation terms. Equation (l) defines the value of Tq^ since in the 

dTg 

limit &e t «o, rr“ 0 equation (l) reduces to 

dt 

|jDg + - Tg « 0 (10) 

For any specified conditions the term in the hraokets is a constant 
which, once determined, is one of the two constants in an equation of 
the type aX^ + X - h = 0 * The positive real root of equation ( 10 ) 
may generally he located qiilckly hy trial and error or hy Newton's 
method of approximation as given in a numiber of standard mathematical 
works. 


KESUITS AND DISCDBSION 


Validity of skin-temperature calctOatlons . - Skin temperatures 
have heen measured at two points on the conical warhead of V -2 mlsalle 
numibar 21 (reference 8) . These flight measurements can he used as an 
index of the applioahility of equation (3) to flight and of the kind 
of errors to he expected when equations (l) and (3) sure used for skin- 
temperature cel-culations under flight conditions. 

Based on data supplied l:y the Naval Besearch Lahoratory, the values 
of skin temperattire measured on hoth the fcarward a1.um1.nuTn-alloy section 
and the rear steel section of the warhead of missile number 21 are 
given in table lU along with the radar-tracking data, the pressure 
and temperature of the amhlent atmosphere, and the pertlrasnt structural 
details. The values of pressure were determined from Instruments 
located in the missile; the amiblent temperature tp to 50,000 feet was 
determined from radiosonde data, and temperatures above 50,000 feet 
are as reported in reference 8. 
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The tiiDB history of slcln tenperatvcre for the alumlmaa-alloy and 
the steel sections has heen calculated hy the niethods outlined in the 
present report. For the aluminum section, a raJLue of skin factor of 0.34 
vas used, corresponding to a thickness of 0.105 inch. This -mlue 
vas selected since the dimensions of the temperature gage attached to 
the Inside of the skin were so large in comparison vlth the skin thickness 
that it vas judged that the ten^ratures measured vere more nearly 
representative of the increased thickness at the gage station. Secause 
of the much greater heat capacity of the steel section, no correction 
vas made for the loceLL Increase at the gage station, hut the variation 
vlth tenqperature of the akin factor for steel, as shown in figure 2, 
vas allowed for in making the calculations. To allow for radiation, 
a value of e » 0.9 vas used as representative of the emlssivity of 
the laoq.uer finish of the nose section. Since the effects of solar 
heating vere found to he negligihle, no atte]]^>t vas made to allow for 
them. 


The measured skin temperatures for this 7-2 missile are shown 
in figure 3(&) • Some of the f 11 ^t -path data are shown in figure 3(h) 
vhich, for purposes of conparlson, also shows the skin temperatures 
measured on the aluminum section. Also shown in figure 3(a) are two 
calculated time histories of skin tesperature for hoth the aluminum- 
alloy and the steel sections. For the solid line in each case the 
value of h vas based on 1^, -for the dotted line the value of h 
vas based on Tg. 

The preclslm of the measured skin temperatures Is given as ti8°p. 
(See table IH.) Considering both the aluminum-alloy and liie steal 
sections, the agreement between measured and calculated skin temperatux^s 
is better for the calculations in which h is based on Tg. The 
measured temperatures lie between the calculated temperatures in the 
case of the aluminum-alloy section but are greater than either of the 
sets of calculated temperatures in the case of the steal section. 

In order to make the calculations shown in figure 3(a), eqtxiation (3) 
has been used whether or not the value of Eq based on Tg fell in ^e 
range of the experiments of reference 2, that is, 2 X 10^ to 2 X 10°. 
Although at 100 seconds the value of Be for the aluminum section had 

fallen to 106 f dvring that portion of the flight in which significant 
aerodynamic heating effects oooiirred, 10 to 70 seconds, the veLLue of Eq 
for both the aluminum and steel sections fell in the ranfp lo"^ to 
2 X 10^, values which do not differ from the range of the experiments 
of reference 2 by more than a factor of 10. In view of the agreement 
obtained it is felt that the following remarks are Justified. 

1. The experimental relation for heat -transfer coefficients for 
conical bodies under supersonic conditions, 

Nu = (0.0071 + 0.0154P^*^) Bq®*® 

glTOii “by E“ber In reference 2 sisy "be extended to flight conditions . 
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2. Satisfactory agreement "between measui^cL skin temperatures and 
temperatures calculated on the "basis of Eher’s equation h as "been 
o"btained in the case of one fli^t o-ver a Mach numher range up to k-,6 
and altitudes up to 250,000 foot. 

3* In the extension of Bber’s equation to fli^t conditions 
"better agreement "between measured and calculated skin teitqperatures in 
the case of Y-2 missile numher 21 was obtained when the "boiuadaiy- 
layer temperature was used as the characteristic temperature than was 
obtained when the skin temperatince was used, but one experiment is 
probably not enottgh to make a final conclusion about the best choice 
of temperature. 

i)-. In the absence of further experimental data the characteristic 
length to "be used in Eber’s relation may "be teken as the distance of 
the point under consideration from the vertex of the conical body, 
measured along the stxrface of the cons . 

Bangs of skin temperature . - In view of the agreement "between 
msasvtred skin temperatures on "Y-2 missile number 21 and the temperatures 
calculated by use of equation (2) and equation (3) "with h based 
on 1^, a series of general calculations has been made. These calculations 
are designed to illustrate the rangs of skin temperatures which can 
be encountered in the design cxr operation of supersonic airborne conical 
structures and to llltistrate the relative significance of the six "basic 
parameters in determining skin tenperature. The following values of 
the basic parameters wei?e selected for investigation. 


E 

( ft ) 

M 

^ Btu > 

l 

( ft ) 

€ 

P 

(deg) 


0 

1 

0.05 

0.1 

0 

mm 

35,000 

2 

.2 

.5 

0.05 


80,000 

3 

.5 

1.0 

.2 


150,000 

5 

1.0 

2.0 

.5 


250,000 

8 

1.5 

5-0 

1.0 

Mam 


Time histories of skin temperature and equild"hrium skin temperatures 
were computed. The time -history calculations are for the convenient 
eissuired flight, plan in which the altitude remains constant, the value 
of M is reached instauitaneously at t = 0, and the initial skin 
temperature is eqtial to In order to simplify the presentation and 

to redtioe the number of calculations, the intesmisdlate underlined values 
of the basic parameters were selected as representative. Each parameter 
was investigated for the five values given, while the other five parameters 
were held at the underlined representative values. 
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The resxilts of the time-history calciLLatioiia are given in figure k, 
vhlch shows "ttie change ih the time history of skin temperature for the 
various values of the hasio parameters . The solid linee^ marked day 
conditions, include the maximum effect of solar heating. For 
dotted lines, representing night conditions, the value of C was 
taken as zero in evaluating equation (2) . 

Since the time of flight of most superscmic missiles and airplanes 
is likely to he brief, figure h shows that the skin tearperatu 3 ?e rise 
during a limited time of flight, with the proper selection of basic 
parameters, may be held to structurally permissible values. In 
particular, the value of Ts which will be reached in any specified 
time at any specified Mach number is smaller at high altitudes, smaller 
for large values of Fg and I, and smaller for low values of p. 

For night conditions, increasing the value of € will always ies\ilt 
in lower values of Tg. For day conditions, however, solar radiation 
can result in more rapid heating so long as the value of Ts is leas 
than about 700® F abs. The effect is small^ however , and after the skin 
teu^erature has risen hi^ enough so that the net radiation exchange 
results in a cooling of the surface, an appreciable reduction of skin 
temperatures can be achieved by employing a surface emlselvity that 
approaches a value of 1.0. (See fig. 4(d).) 

Since the rate of increaee of Tg may become very small as Tg 
approaches its final or equilibrium valtie Tq, the time histories of 
skin temperature shown in figure 4 have not been extended beyond the 
time at which the initial ten^rature difference has fallen to about 
10 percent of its initial value. The value of Tq is shown, however, 
in figure 5 for the sane combinations of the basic parameters ae were 
used in figure 4. The values for both day and night conditions Illustrate 
the maximum effect of solar heating, which is seen to be small in all 
cases except at veiy high altitudes. For conrparison, the value of 
la also shown in figure 5- Values of 0^ were confuted using equations (6) 
and (9) and reference (3)) values of Tq were computed using eqmtion (10). 

Skin temperatures for a typical supersonic airplane .- In order 
to illiistrate the magcitude of ajerodynamic heating effects for 
practicable flight conditions, the time history of skin temperature has 
been computed for three representative flight plans for a typical 
superscaiio airplane. The results, shown in figure 6, are for a point 
1 foot back from the apex on a conical nose of total included angle 0 

equal to 30®. Eepresentative values of € =0.2 and Fe = 0.2 

f t2op 

have been chosen. lEhis latter ved.ue is applicable to 0.034-lnch- 
thlok stainless steel, or 0.06l-inch-thiok 24S-T. 
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Flight plaa A (fig. 6(a)) is a hlg^-altitude run, foUowsd "by a 
constant” speed, gl l de to lower altitudes . The skin teiiiperat\ire does 
not exceed 5^^6° F ahs. (86° F) . Flight plan B (fig. 6(h)) is a run of 
2 minutes duration at hi^ speed and constant altitude . Because of the 
long time spent at high speeds, the skin temperature rises to 8l|-7® F ahs. 
( 387 ° F) . The peak in temperature closely follows the peak in speed. 
Flight plan C (fig. 6(c)) is a short hurst of speed to mw-xinnijffi Mach 
number, with full power, and a longitudinal acceleration of approxi- 
mately 2.8g. Because of the high speeds involved, the tenqwrature 
rise is very rapid, and continues for approximately 20 seconds after 
the power is cut off. At this time the value of T 3 has fallen helow 
the value of Ig . The maximum teniperature reached is 830 ° F ahs . 

( 370 ° F) . Althovigh the skin temperatures reached in the high-altitude 
run, plan A, are not excessive from a structural standpoint, the skin 
tenqperatttres reached in plans B and C are excessive. Calculations 
for plan C indicate that an Increase in skin factor from 0.2 to 0.5 
would result in a decrease in msulmum skin temperature of approocl- 
mately l40° F to a value of 69 O® F ahs. (230° F) . 


COHCLOSIOMS 


1. The experimental relation for lieat -transfer coefflolonts for 
conical bodies imder supersonic conditions 


ifu = ( 0.0071 + 0 . 0154 p 0 . 5 )Bq 0‘8 


given hy Eher in reference 2 may he extended to flight conditions. 

2. Satisfactory agreement between measured skin temperatures 
and temperatures calculated on the basis of Fber's equation has been 
obtained in the case of one flight over a Mach number range up to 4.8 
and altitudes up to 250,000 feet. 

3 . In the extension of Eber's equation to flight condition^ better 
agreement between measured and calculated slcln temperatures in the 
case of V-2 missile number 21 was obtained when -the boxaidary -layer 
tenporaturo was used ae the characteristic tenperature than was 
obtained when the skin tenperature was used, but one oaperiment is 
probably not enou^ to make a final conclusion about the best choice 

of temperat^^re . 

4. In the absence of further experimental data the characteristic 
length to be used in Eber’s relation may be taken as the distance of 
the point under consideration from the vertex of the conical body, 
measured along the surface of the cone . 
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5* The skin tempera tvires which will he reached on supersonic 
missiles or airplanes may he very much less than either the eq^iilllhrium 
skin temperature or the boundary -layer temperature, at hi^ altitudes, 
and eJ.so at low altitudes if the fll^t duration is short." 

6. By proper selection of the basic parameters, the skln-temperatvire 
rise during a limited time of flight may be held to structurally 
permissible values. 

T* Low values of skin temperature, and low rates of increase of 
skin tempera tui*e at any given flight velocity are associated with 

(a) High altitude 

(b) High val-uea for the skin of the product - specific heat x 

thickness x specific weight 

(o) Hl£^ values of emlsslvlty 

(d) Small cons angle 

(e) Largsr distances back of the nose. 

8. The effect of solar heating on the skin teii 5 )erature of supersonic 
airborne structures is small at altitudes below 150,000 feet. 

9. The radiation received from space and the outer atmosphere 
other than solar radiation, may usmlly be neglected in the calculation 
of skin temperatures for supersonic eiirbome structures. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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AEEEBDIX A 


DEEIVACTON OF !EHB EIFFEREIIPriAL EQUATION 
FOR SKIN TEMPEEIAIPURE 


The qiuantlty of heat transferred In unit tin© throu^ a fluid 
film Is glTen as 

dQ » m(^ - . (Al) 

The heat lost hy radiation Is given hy 

dQ = A«oTg dt (A2) 

and the radiant heat ahsarhed from space and the outer atmosphere Is 
given hy 


dQ 


A€ct8Ta dt 


(A3) 


Ec[uatlon (A3) Is hased on measurements mads on the ground at different 
angles to. the zenith reported In reference 9 where It Is shown that the 
radiant energy received from, space and the aiaiosphere can he correlated 
vith the local ambient atmospheric temperature if the atmosphere is tahen 
as a black body of temperature and emissivi-ty 8. The quantity 8 

varies linearly with the square root of the proportion of the total 
atmosphere included In the path of the measuroments . Since the proportion 
of the total atmosphere above any point to the proportion above a 
point at sea level is given by the ratio p/Po, ^ has been taken as 

proportional to /p/po in ihe present report. Although this extension 
to high altitudes of measurements made on the ground Is admittedly an 
approximation, the validity of which must await experimental conflimatlon, 
the approximation Is conflldered a better one than would be the case 
If some nominal black -body ' temper ature were adopted for 'the a-fcmosphere . 

The variation of 8 with /p/Po Is shown in figure 7, as adapted 
from reference 9- 


The equation 
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can l3d used to represent the maximum amount of radiant heat ahsorhed 
from the sun when the svirface of the skin is perpendicular to rays 
of the sun. The value of the solar constant C is given In reference 10 
as 1.94 cal/cm2 min (0.1192 Btu/ft^ sec). Since the solar constant is 
defined as the average total energy received from the sun per unit area 
per unit time at the top of the earth’s atmosphere, corrected to the 
sun's mean distance from the center of the earth, eqtuatlon (Ah) 
neglects the effects of season and the variation of atmospheric 
ahsorption with altitude . The emlsaivity € is also used in oq,uatlon (A4) 
Instead of the ahsorptlvity for solar radiation. While the ahsorptlvity 
of some materials of aircraft construction is hlgier than their 
emissivity, equation (All-) is felt to he a valid approximation for the 
maximum effect of solar radiation particularly for the reuigo 0.2 < e <1.0, 
A more complete consideration of the effects of solar radiation, 
particularly as to the Influence of selective ahsorption and emission 
on surface tenqperature, is heyond the scope of the present report. A 
brief discussion will he found in reference 11. 

The heat required to raise the temperature of a hody is 


dQ = cV dT (A5) 


It is assumed that (a) the heat transferred to the interior of 
the cone can ho neglected since the skin is hacked hy air only, (h) 
that radiation losses to other parts of the structure are negligihlo, 

(c) that radiation from the inner surface is halanced hy radiation from 
the opposite side of the cone, (d) that the thermal coxiduction of heat 
away from the point under consideration due to a tenqperature graulient 
edong the surface is negllgihle, and (a) that the skin is sufficiently 
thin that no temperature gradient exists in the skin perpendicular to 
the surface , duo to the flow of heat at the external surface . Under 
these assun^itlons, for tlie portion of a cone perpendicular to the solar 
radiation, the following equation, comhlnlng equations (A5), (Al), 

(A3), (A4), and (A2), would hold 

^ ^ “ Ts) dt + AeoSTj^^ dt + AeO dt - AaoTg^ dt (A6) 


Since the weight of the section of skin of area A can he expressed as 


W = Atw 


(AT) 
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IT 


and since In eguatlon (a 6) dT * diTs , the following hasic differential 
eqmtion for skin ton^^ratiire is obtained 

dTg h(TB - T|^ + 
dF" “ oTw 


This eq.viatlon is given in a convenient general form, in order to 
permit a definite evaluation of the contribution to dTg/dt of the 
convection, reuilatlon, and solar terms. In cases of practical Interest, 
a brief studjr will show which of the terms in eq^uatlon (l) may be 
omitted without the introditction of significant errors. 
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APPENDIX B 

NUMERICAL CALCULATIONS WITH THE HEAT BALANCE E»UATIQN 


Niimerlc^ Integration . - In tlxe stop-by-step numerical Integration 
of eqmtlon (l) it has been fotuid convenient to rearrange the terms as 
follows: 


^ - BTs^ - Ts 

dt 9 


(Bl) 


where 


B = ¥ (°F abs.)"^ 
n 



The magnitude of B is a meeisuro of the relative significance of 
the radiation and conduction terms. The quantity Tp which may be 

termed the potential temperature is a fictitious ten^erature resulting 
from all of the heating terms in equation (l) . The quantity 9 is 
the time constant of the thexmeO. system and is a measure of the time 
required for an initial ten^ratiire difference to fall to l/e times its 
original value (when radiation may be ignored) . 

Numerical Integration over successive time Interveils t]_, t 2 • • * tjj, 
^+1 when neither the speed nor altitiide is changing is a simple 

iterative arithmetic process which can bo qtiiclcLy carried out. For 
fli^t paths in which altitude or speed vary with time, the process is 
essentially no more complicated once the average values of Tg^ h 
ani^ 6 Ta,^ have been determined for each time interval. For such * 
varying flight paths, equation (3) would be written 



(B3) 
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whare 



Althou^ the temperattire Ts^ which, will ho reached at tha end of 
the time intenral At is not known in a step-hy-etep integration, a 
reasonahle value may he assumed. The valiae of Tm corresponding to 
this assumption is then tised to calculate a value and this first 

calculated value is con^jeirad with the assumed value. If there is a 
diffeirenoe between the two values, rapid convergence will generall; 
obtained if a new assumed Tsg is used which lies midway between the 
assumed value and the first calovilated value . Two or three iterations 
of this type will generally give assumed and computed values of Tg^ 

which agree within 0.2° F, a value which is judged to be sxafficlently 
close to prevent serious accumulative errors. 

Approximate integral solution .- Equation (Bl) can be Integrated 
approximately by analogy with tha integral of the basic equation for 
heat flow, generally known as Hewbon’s law. The application of this 
law to thermal problems and. In particular, to thermometers has been 
given by a number of authors. (See, for example, reference 12.) The 
tenperature of a thin sheet of skin initially at tenperatvtre Ti which 
is plungjsd into a bath of temperature T^, from which it is insulated 
by a layer with heat -transfer coefficient h, is given as a function 
of time by the equation 
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The solution of eq^uatlon (b 4) is 


% - T 
- Tl 


_t 

= 00 


(B5) 


If the term Tp - BT^ in equation (Bl) can he oonsldeired a constant 

over a small time interval At, then an integral of equation (Bl) 
can he -written as 

At 

= o "ST (b6) 


where, as before, Tm is -the mean value of T over the time interval At. 
Equation (b 6) may he re-wrlt-ten in a farm, siiltahle for s-fcep-hy-s-fcop 
integration as 



(Tp - BTm3 (l - 9‘e^ + Tie" 0 


Since Tm over -the time Interval At is not loiovni, Ti can ha used 
as a first assumption for Tj^, and the resul-tant value of T 2 used 
-to make a het-ber assuo^tlon. In view of -the assumptions used in the 
derl-vatlon of equation (B?)# -valxies of At/9 should not exceed 0.25* 

A sample calculation using equation (b6) Is given in -bahle IV. For 
a rapid approximation, the 1-fceratl-ve process may he eliml.na-ted and -the 
value of T 2 obtained by assuming that T^ = Ti may he used as the 
■value of- Ti for -the next time In-terval. 
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TABLS I 


COMPARISON OF STACHIATION TEMPERATOEES CCMPDTBD 


BY TJSHSi CONSTANT AND VARIABLE SPECIFIC HEAT OF AIR 



^A 

a 392 *^ 

^A 

*= 518.4 

H 

^ . 
7 = 1.4 

eg.uation (9) 

^ . 
7 = 1.4 

T^ 

equation (9) 

1 

470 

470 

622 

622 

2 

706 

705 

933 

929 

3 

1098 

1087 

1451 

1418 

4 

1647 

1600 

2177 

2060 

5 

2352 

2210 

3110 

2844 

6 

3142 

2860 

4251 

3770 

7 

4232 

3760 

5599 

4830 

8 

5412 

4680 

7154 

6 o 4 o 


Temperatures In °F ats. 
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TABLE II 

THERMAL EROPHRTIES OF MATERIALS 


Specific heat and apeclfio weigjit^ 

Material 

O 

(Btu/lb °F) 

(ll>/ft 3 ) 

Aluminum 



Piare 

0 . 226 ( 0 - 100 ° C) 

168.5 

2te“T 

.226(0-100° 0 ) 

172^2 

WroTig^t iron 

.1138 

483.6 

Stainless steel 

.142 

492.5 

Magnesium 

- — 

. 249 ( 0 - 100 ° C) 

108.6 

Specific heat of iron® 

T 

« T 

c 

(°c) 

(°P abs.) 

(Btu/lb °F) 

0 

492 

0.1055 

200 

852 

.1282 

400 

1212 

.1509 

6oo 

1572 

.1737 


^Data from reference 6 . 

*Data taloen from table entitled ^Specific Heat - 
Variation with Tea5>eratxire “ of roferenoo 7 * 
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TABLE III 


SKER-TEMPlEATUEffl MEASUREWEHTS- 
V-2 MISSHE HIKBER 21 FLRSD MARCH 7, 19^<-7 
MTA SUPPLIED BY NAVAL RESEARCH LABORATOBY 


Time after 
launching 
(sec) 

Altitxide 

(km) 

Velocily 

(m/soc) 

P 

(mm Hg.) 

(°K) 

T 

Aluminum 

(°c) 

^Steel 

(°c) 

0 

1.2 

0 

61^-5 

280 

22 

30 

10 

1.7 

109 

600 

279 

22 

30 

20 

3.5 

260 

480 

267 

22 

30 

25 

5.0 

348 

410 

261 

22 

30 

30 

7.0 

438 

300 

246 

25 

30 

35 

9.4 

530 

210 

228 

35 

31 

40 • 

12.3 

646 

130 

218 

45 

40 

45 

15.9 

810 

76 

215 

62 

55 

47.5 

18.0 

900 

52 

215 

77 

65 

50 

20.4 

995 

33 

215 

87 

75 

52.5 

23-0 

1090 

23 

215 

97 

84 

55 

25.9 

1200 

14.5 

216 

108 

92 

57.5 

29.0 

1305 

9.6 

216 

118 

98 

60 

32.4 

1420 

5.8 

217 

124 

105 

65 

: 39.9 

1540 

2.3 

270 

135 

116 

70 

: 47.5 

1490 

.95 

315 

142 

i 122 

75 

54.8 

1440 

.40 

320 

147 

122 

80 

' 61.8 

1390 

.18 

299 

I 156 

125 

100 

87.5 

1200 

,003 

a)0 

! 160 



t.5 




±10 

±10 


Cone angle, P, 26°. 

TMcknesB of 3S ed-umlniam -alloy section., O .091 Inch. 

Thiolmess of Spec. 4055, grade M steel section, 0.109 inch. 

Size of teisperature gage, 0.014 inch X 0.6 inch x 3*5 inch approximately. 
Location of gage on alvnnimcn -alloy section, I = I.5 foot. 

Location of gage on steel section, I = 2.6 feet, 

Svirfaoe finish, orange laoq.uer. 




M p 3 

H . X50,ooa n 
p « 30° 
l o 1 ft 

* - 0.2 

Pg - 0.2 Bttt/ftF °P 

■ 373»5° ^ ®*®* “* t a 0 


nCBlJC IT 

aAMPI£ TW’SISECailT QALBOLUnaa 


1441° P "B8. 

0 > 800 

3.8 X 10"* Btii/rt2»o ftliB. 

At ^ 60 «oo 

2.904 X 10'^ ^ Blia.'^ 

^ - 0.U4O25 

10.02 X 10^ “P aDfl.*^ 

A- 


a 0 - 0.09224 

1531" P 0l)B. 



1 - a ® = 0.10776 


a 0.10776 + 0.8922 >ITi 


(3) f (10) 

& 


573.? 573.5 
573.5 


I MV ”p " ^ 

(5) X B 1531 - (6) 0.09224 X (3) O.IO776 x (7) (A) + (9) 


0.1003 X loi^ 
.15122 
.IJUfl 


1503.9 

1493.1 

1493-2 

l4T9.fi 

1465-2 

1463.7 


1075.6 

1062.1 

IOGL.6 


1.3304 

1.3711 

1.3^ 


3J-95.9 

US7.7 

H66.3 


IO6B.6 

1007.7 

1087.7 


660 1007.7 
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M 

Figure /.-Temperature recovery factor 

as a functfon of Mach number and cone 
any/te . Data from reference ^ . 
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^^abs 



(a). Ca/cuJahd and measured skin iemperahre • 

Figure s , - Time hishry of calculated and measured ej^uantities 
during ascending portion of flight path \F^ missde No. Zt . 
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sec T 1^ ^ 7 ^^. Ve/ocit^y 
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T~ime , sec 

(a) Mach number M . 

Ffgfure Change /n skin - femperai^ure iime h/siory 
for various va/ues of desfyn parameter . 




Skin femperaiure 
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£ to 15 20 

Time , sec nt/n 

M Mfifude , M, 

Figure ^ -Conh'nuea/ . 


Skin femperaiure , 



3kin femperai^ure 
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Figure 4 - . -Coniinueoi . 


Skin temperafure 
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14 ? 


Figure 4,- Coniinued . 


Skin temperature 
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(e) Length , I . 






Figure 4 Continued . 


Skin temperature 
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(f) Cone angle ^ ^ . 


F/gure 4 - ~ Concluded . 


Skin temperature 
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H= so . 000 _L 

0,Z 


Z - /f^ _ 


^=30"* 


Day 

Ts 

Night 


~^~3 r 


(a) yanafion i/^^irh Mach number, 
rtgure 5.- l/arfat/on of equthbr/um skin fenop 
aiure Jg ; c/ay and nig/^/ conc/r/iona , 
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Figure J . - Conh'nued . 
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ic)Varfcit/on tNifh emnis^ivii’y. 





(d) Von a ’t ion uw/'/’h lengi-h. 



<?jdeg 

(.e) Vanahon k\nfh cone angle. 




Figure 5 - Concluded^ 





n^ure e.^T/me hisfoty of Maeh number ^ Mwh ,und oahuhfed skin hmperahre for 
fHghi plans of a fj/pica! supersonic airplane . }” t f! j Fs*0,sBluJst^ ff *F , » 
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O / _2 3 4 3 6 

Time j mm . 

(b) Fhohi plan B . 

Figure O. - Continued. 
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Mach number 


( 


— / 60 \Xl<f 


4 : 

^ >-900 




n Mach number — ^ 


Q) Q) 

5 :s 

'^90 TOO 

$ ^ 

a 

h 5 «w 


S soo 


A 


'^/IlhiudeM 



^mperaiure 


Figure - Concluded, 


& a 

Time ^ mm 
(P) Flight plan C. 
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«♦ 



Figure 7 ".- Sky radiation factor . Daia from 

reference 3 . 
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